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ABSTRACT 

This paper brings together a number of algorithms for the estimation of the envelope and 
instantaneous frequency functions suitable for machinery diagnostics and structural health 
monitoring. This includes the classical approach based on Rice's frequency, the most 
commonly used method based on the Hiibert transform and recent developments related to 
the Wigner-Ville distribution and wavelet analysis. The study includes two application 
examples from fault detection in gearboxes and identification systems, The paper shows that 
the envelope and instantaneous frequency functions are an important link between classical 
Fourier approach based methods and time-variant analysis. 



1 A 
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1. INTRODUCTION 



The viewpoint of this paper is that there exist three major areas within machinery 
diagnostics and structural health monitoring where the instantenous characteristics, namefy 
the envelope and instantaneous frequency, are used. These are: damage detection in rotating 
machinery, identification of nonlinear systems and deterministic approach to stationarity. 

Machines and structures in operation generate forces and motions that produce 
different types of vibration. Any fault in machine operation or structural damage can be 
considered as an additfonal excitation which results in a change of vibration response. Often 
the overall process of vibration response exhibit modulations, Intermediations or nonlinear 
behaviour due to operational nature or fault. Identification of modulation sources can provide 
valuable information of many faults In rotating machinery, particularly in gearboxes and ball- 
bearings. Modulations In gearboxes are mainly caused by the tooth meshing process, due to 
varying stiffness, or by the eccentricity and local tooth faults 1 . This is exhibited by modulation 
sidebands in vibration spectra, 

Intermodulations results from sum and difference frequencies of low frequency 
harmonics of the shaft speed, harmonics of the tooth meshing frequency and modulation 
sidebands. Modulations in ball-bearings are mainly caused by local surface defects which 
produce impulses often leading to resonances. For advanced defects, simple parameters 
such as kurtosis or crest factor 2,3 are sufficient to detect damage. Often modulation and 
intermodulations which result from damage can be identified by means of ZOOWI-spectrum 
and cepstrum 4,5 . For complex machines the signal demodulation procedure, or in other words, 
the analysis of the envelope and instantaneous frequency functions. Is very effective to 
identify local tooth damage 6 ' 113 and local defects in ball-bearings 11 " 14,9-10 . The amplitude 
demodulation procedure is often called envelope detection or, when used in ball-bearings, the 
high-frequency resonance technique. 

The second important area where the envelope and instantaneous frequency 
functions are often used is identification of nonlinear systems. Since nonllnearities may result 
from damage in a system, the distinction between linear and nonlinear behaviour is an 
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important problem in damage identification. Many different procedures for the identification of 
nonlinear systems from vibration test data have been developed in structural dynamics. It is 
well known that many types of nonlinearities cause a varying nature for the restoring forces 
and natural frequencies of the system. This varying nature of vibration can be studied using 
the envelope and instantaneous frequency functions which lead to nonlinear back-bone and 
damping curves. 

The overall vibration produced by machinery and structures is either stationary or 
nonstationary. The analysis of the envelope and instantaneous frequency functions often 
forms the deterministic test for stationarity In vibration analysis 15 . In this approach, the 
process is said to be stationary if it consists of components such that their instantaneous 
behaviour, as described by the envelope and instantaneous frequency, does not depend on 
time. Conversely, the process is non-stationary if its instantaneous characteristics depend on 
time. The instantaneous characteristics are used in machinery diagnostics to detect local 
faults. This fault detection method is in fact a test of stationarity and any departure of 
instantaneous characteristics from time independent functions is considered as a symptom of 
a fault in the system* 

Recent developments in signal processing show that the envelope and instantaneous 
frequency functions are an important link between classical Fourier approach based methods 
and time-variant analysis, as shown in 15 . There exist a number of procedures, developed 
over the last twenty years, which can be used to estimate the envelope and instantaneous 
frequency functions. It appears that the Hilbert transform based procedure 16 is the most 
common approach to the problem. Recent developments include the Wigner-WIe 
distribution 17 and the wavelet transform 1 *. However, there exist a number of different 
aigorithms which can be used for single- and multi-component signals (frequencies), 

The arm of this paper is to bring together the algorithms of signal demodulation. It is 
hoped that this form of presentation will help to implement the most suitable approach for a 
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given vibration problem. However, the paper does not intend to survey various application 
examples. 

The structure of "the paper Is as follows. Section 2 shows the complexity of 
modulation and intermodu!ation processes which can be found in rotating machinery. The 
algorithms of the amplitude and frequency/phase demodulation are given in Sections 3 and 4, 
respectively. This is followed by Section 5 with examples related to fault detection in 
gearboxes and identification of nonlinear systems. Finally, the paper is concluded in Section 

2, MODULATION AND I NTERMODU LATION PROCESSES IN MACHINE DIAGNOSTICS 

In machinery, modulation is the variation in the value of some parameter which 
characterizes a periodic osciliation. This process is an effect of either design features of 
machines (e.g. crankshafts) or more often a result of a fault Thus the identification of 
modulation sources can provide valuable information about possible faults in machinery. 
Amplitude modulation processes, in general, are caused by impacts and Impulses in the 
structure, generated by different phenomena, for example, cracks or missing teeth in gears, 
point defects in ball-bearings, blade cracks in fans or any other local faults. The impulses are 
repeated periodically for each revolution of the wheel, shaft, etc. The frequency modulation 
processes are an effect of parametrical excitation in the structure (e.g. fluctuation in tooth 
loading in gearboxes, varying stiffness of teeth in gears, varying cross-section of shafts, etc). 
More details can be found in t,19 ~ 25 . in what follows some spectral properties of modulated 

processes will be relevant to machinery vibration as summed. 

> 

The vibration signal generated by a machine can be analytically described as a compiex wave 

*W=Z&+*r WV^ ^ {t)] + *(') 
r=l 

(1) 
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where a r (t) is the slow-varying amplitude modulation process, £.(r) is the fast varying 
phase modulation process, f or are the carrier frequencies (o 0 r=2xf or = const) and n(t) 
is the noise. Amplitude and phase modulation processes can be represented by, 



k=\ 



(2) 



TY). T 



(3) 



where M ir are the coefficients of amplitude modulation intensity and, J3& are phase 
deviations. Substituting (2) and (3) into Eq. (1), it follows that 



(4) 



Eq. (2) can be replaced by 5 



(5) 



The welt known Fourier-series expansion, 



(S) 



where J n (Pfr) is the Besse! function of the first kind, together with Eqs, (4) and (5) yields, 



&m\$£tAxJ | Q< _ [)[q 2 0 0 PM 4 ' 3 5 ft * s ^ ,fto * b ''HARRISON GODDARO F qq^^-^w* ^ u *^^"<^^ ^ q q . 



2« 



+ 0O 



(7) 



J„ ) in Eq. (6) is the Basset function of the first kind. Eq. (7) can be written as, 



(3) 



or finally, 



n{t) 



(9) 



Analysing Eq. (9) one can notice that in the frequency domain the vibration signal 
x{t) is represented by R families of spectral components. Each family yields five types of 
side-frequency terms: 

A i €^ orl ~ carrier frequency 
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A 2 e J ^ or±Qir ^ - side frequencies due to amplitude modulating waves 

A 3 e X% ±ff %)' - side frequencies due to phase modulating wave 

A 4 e J ^ a} ° r±nmiy±n<x>2r ± ±na} & ) f - beat side frequencies due to phase modulation 

As e j{^ r ±a lr ±n 2r ±^ ±& ir ±n<^ r ±nc^ r ±„±n<0} it )t _ beat sjde f reqaen cies due to 

amplitude and phase modulation, 

where A t are the Bessel function components according to Eq. (9). Fig. 1 shows an example 
of the amplitude and phase modulated signal, which includes two families of spectral 
components, and its spectrum. 

In the case of the pure amplitude modulation process =0)the total energy of the 
signal will increase. The energy of the carrier remains unchanged and the additional energy 
which comes from the modulation wave will be represented in sidebands. This process yields 
only one pair of sidebands: sum and difference components. In comparison the total energy 
of the phase (frequency) modulated wave remains constant To achieve this the energy of the 
carrier will decrease and the difference will be distributed to its infinite number of sidebands. 
In addition phase (frequency) modulation can produce beat components, which are not 
characteristic for the amplitude modulation. In practice, for example, the carrier frequency in 
the gearbox deals with the meshing frequency. Thus one should remember that greater 
advancement of a fault does not need to cause the increase of the level of the meshing 
frequency. This means that the level of the meshing frequency might not be a good indication 
of the advancement of a fault. Additionally, if one considers the phase angle between the 

amplitude and phase modulation waves i.e. a r (t)=^M ir CQs(Q. i7 .t+Q. iro ) , thus the 

asymmetric distribution of sidebands will be observed. These brief considerations show that 
the spectrum of the vibration signal generated by a complex machine might be very 
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complicated. Fig, 1 shows an example of amplitude and phase modulated signal and its 
spectrum. 



3. AMPLITUDE DEMODULATION ALGORITHMS 



Amplitude demodulation is the process of extracting slow variation of the amplitude 
from the modulated signal. This operation falls into two broad categories of synchronous 
detection and envelope detection 26 . Both categories are summarised below. Additionally, the 
algorithm involving Discrete Fourier Analysis (DFA) 27 is also described. 



3.1. Synchronous Detection 



Synchronous detection uses the operation of frequency conversion. The concept is 
outlined in Fig. 2, Here the oscillator is assumed to be exactly synchronised with the carrier 
frequency. 

Assume an amplitude modulated $\gna\,x(t)=a(t)cosco 0 t where co 0 is the carrier 
frequency. Multiplying both sides by 2cos<o 0 t one obtains 



2x(t)oQ§a? 0 t = a(/)cos 2co 0 1 

(10) 



The first term on the right side is the one we want, since it is an amplitude modulation 
signal. The second term can be removed by low-pass filtering. This procedure requires the 
carrier frequency determination in order to translate the signal spectrum to a new carrier 
frequency equal to zero. Nevertheless some degree of asynchronism must be expected, 
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since the carrier frequency cannot always be exactly determined (for example in the case of 
varying rotational speed of the gearbox), This means that instead of multiplying by 2cos^ t 
one may use the value 2cos(^ which results in, 

2*(*)cos(<» 0 +G>' 0 )l " a{t)to$<o' 0 f+^(?)oos(2^ +0 e o )t 

(11) 

Here the second term can still be removed but the harmonic relationships in the 
modulation signal aif) are destroyed. Therefore Instead of a(t) one receives a(t)cos&' 0 t , 
which in the frequency domain, according to the modulation theorem 28 , gives us the shifted 

spectrum without any possibility of correction. The second difficulty with this 

procedure is related to the digital filtering operation which reduces the efficiency of the 
algorithm, 

An alternative implementation of this method is presented in section 4.2, where the 
frequency modulation signal tp{t) is assumed not to be equal to zero, which allows one to 

compute either a(t) or #(t). 

3.2. Envelope Detection 

There exist a number of envelope definitions. All of them come from random vibration 
or communication systems theory. Physically, for the narrow-banded process x(t), the 
envelope is a smooth curve joining the peaks of jc(r), A few envelope definitions and their 
statistical properties have been described by Langley 29 , In what follows a few possible 
implementations of envelope detection methods are briefly discussed. 

3.2.1. Local maxima distribution 



2001 14:37 " HARRISON GODDARD FOO 



NO. 016 



P. 



The envelope as a local maxima distribution was defined by Crandall 30 " 31 as being a 
smooth gradial curve joining the peaks and the expected time spent by the envelope between 
the level a and a+da Is just the number of those peaks, whose magnitude lies between a 
and a+da multiplied by the expected period for cycles of amplitude a, in practice the 
method leads to running averages of signal local maxima. If the modulated process is written 
as x(r)=<z(/)cos[d)0 t+#(t)\ the running averages can be carried out as 32 



e{t)~ E r {Max \x(t}}=~ J Max{a{l)QO$[&o * + 



T + t 



r 

(12) 

where E T Is a symbol of averaging. If the phase process is constant ($(t)= const) , far the 
optimal averaging time r» — 6 Eq. (12) can be replaced by 



e(t)=- \a{t)dt~a(t) 

T J 



T 



(13) 

This definition implies the bandwidth restriction. The efficiency of the method depends 
on the choice of the optimal averaging time. It also increases if the process is more narrow- 
banded. The local maxima finding procedure can be done by the approximate construction of 
tangents structure. Unfortunately the envelope is then given as a series of finite number 
samples, which makes exact tangent points calculations impossible. Additionally, the 
integration procedure requires a few samples in a period to reach a good accuracy 33 . 

3.2.2. Energy based distribution 



The energy based distribution envelope proposed by Cranda)! 30 ~ 31 is related to 
stochastic processes, if the stochastic equation of motion is assumed to be of the form 30 



9^ 
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x+M(x,x)+G{x) = F(?) 

(14) 

where x is a displacement response process, A is assumed to be an odd function with 
respect to x , £ is a small parameter, F{t) is the white noise excitation, the envelope 
function a{t) of the random process x{t) is given implicitly by 30 

(15) 

where Jg(<")^. The right-hand side of the above equation is the sum of the kinetic 

o 

and potential energy per unit mass, in the linear case one can set J3=2££d a , F(x,x)=x , 
G(x)={£>o x and Eq. (15) can be replaced by 

(16) 

where a> c is some constant frequency. Langley 23 has suggested to use co c as the mean zero 
crossing frequency given by <o c = I— , where m n is the n-th spectral moment of the single 
sided spectrum of x{t) defined as 

90 

0 

(17) 

This definition of the envelope does not Imply any bandwidth restriction. It is easy to 
notice that such an envelope follows the peaks of x(t). since the envelope and the process 
coincide at a peak, where x(t)-0 . The implementation of the method can be realized 
following Eq. (16), The effect of the algorithm is shown in Ftg, 3. This requires signal 



24 



'19. DEC. 200ri4:44 <Mft "~ , HARRISOH GODDARD FO(f * " ~~ NO. 016""'" ~"P. 29 



differentiation, which is not a simple task 33 The differentiation can be avoided by measuring 
acceleration or velocity where possible and integrating to obtain velocity or displacement The 
accuracy still requires sampling at over a few times the highest frequency of interest 33 . 

Langley 29 has shown that the mean rate at which the envelope crosses a given level 
with positive slops depends on the 4-th spectral moment m 4 . He indicated that some random 
processes, such as for instance the response qf a linear system to white noise, have a 
theoretically infinite value of m 4 . This predicts that the envelope crossing rate is infinite. In 
practice signals generated by machines have a finite value of m 4 , if not, it is still possible to 
filter the response spectrum at a frequency which yields a finite value of w 4 . 

3.2.3. Rice's envelope 

The so-called classical definition of the envelope has been given by Rice 34 . Writing 
the modulated signal in the form 

n 

(18) 

and selecting a frequency f m called the "midband frequency" Eq. (18) can be replaced by 

n 

(19) 
where, 

ic -Z>* cos K^« + ^)\ 

n 

(20) 
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IS ^Cn&afan-fm)* + 4n (?)] 



The expression 



(21) 

is termed by Rice the envelope of x(t) referred to frequency f m . It has been shown by 

Dugundji 36 that this envelope is completely independent of the midband frequency fm> The 
explicit calculation of Eq. (21) is rather impossible. Thus the Rice's envelope is presented only 
as a theoretical case. 

3.2 A. Hilbert transform approach 

The modulated signal x(t) can be replaced on the basis of the analytic signal 

x A (t) = x(f)+Jx H (f) 
(22) 

where x H {f) is the Hilbert transform of x(t) being 2 



-28 



(23) 



Thus the envelope suggested by Dugundji 35 is possible 
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(24) 

It has been shown 29 * 35 that this envelope and Rice's envelope are equivalent and the 
proof of equivalence does not depend upon x(t) being Gaussian. It is easy to notice that for 
x(t) being harmonic (x(t)=Acosco 0 1) the energy based envelope (16) and the envelope 
defined by (24) are also equivalent However, the central carrier frequency of Dugundji's 
envelope is independent of the carrier frequency. 



Langiey 29 has assessed when the envelope given by Dugundji (24) will follow the 
peaks of x(t). The mean and variance of x H (t) when x(t) has a specific value is given by* 9 



x H 



r \ 

\ m 2j 



(25) 



var 



(26) 



where q is a parameter which measure the extend to which x(t) is narrow-banded 

2 

q =1 

(27) 

Thus a(t) given by (24) will follow the peaKs of x(t) (x=0) if q is small {x ff =d) 
which requires x(t) to be the narrow banded process. It has been shown that the mean rate 
at which the envelope (24) crosses a given level with positive slope is always finite, which 
does not take place in the energy based envelope. 
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The signal processing implementation of the envelope proposed by Dugundji 3 * can be 
realised according to Eq. (24). The Hiibert transform definition given by Eq. (23) can be 



written in the convolution form 



,2a 



n t 



(28) 



Taking into account the signum function, 



sgn* = 



-1 /<0 
+ 1 i>0 



(29) 



and the Fourier transform 28 , 



+00 



F[x(t)]=X(f)= jx(t)e' 2 ^dt 



(30) 



the following is obtained 



(31) 



-Flsgn/]=-rJ-T 



(32) 



Thus Eqs. (28), (31) and (32) yield 

F[x H (t)]=X T (f)=X{f)(-jsgaf) 

(33) 



28 
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where X T is the signal X{f) with shifted phase by ^ for the negative frequency 

components and ~— for the positive frequency components. It means, that the Hilbert 

transform x H (t) of the vibration signal x(t) can be easily obtained by calculating the 
complex spectrum of the signa I (*(/)) , shifting the phase (x T (f)) according to Eq. (33) 
and employing the inverse Fast Fourier Transform (FFT). The computation procedure can be 
realized aiso in another version. Going back to Eq. (22) and using the Fourier transform (30) 
one can find 



X(f)=X(f)+jX T (f)=X(f)+sgxfX(f)= 

(34) 



0 /<0 

x(f) /-o 

2X(f) f>0 



This equation represents the spectrum of the analytic signal and can be very easy 
computed using the spectrum of the signal x(t) . It allows the computation of the envelope 
function directly from Eq, (22) as a modulus of the analytic signal. 

The method presented in this section seems to be very effective and is often given as 
an example of an effective algorithm of the amplitude demodulation 37 - 21 ' 33,9 . It has been used 
in a number of diagnostics applications 37 ' 21 ' 3S ' 9,32 . 



3.2.5. Discrete Fourier analysis (DFA) algorithm 



This algorithm has been formulated and tested by Hsueh and Bielawa 27 . The 
assumption is that the amplitude and phase are very slowly relative to the carrier frequency, 
which means that signals are narrow-banded. If they are additionally periodic, they can be 
represented by a time-variable Fourier series. Thus a typical amplitude modulated signal can 
be written as follows 27 



9Q 
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x(t)=x Q (t)+x mc (f)cos a>Qt+x m (t)s\u&Q i+higherharmonic terms + noise 

(35) 

where tf> 0 is carrier frequency, are th€ time-variable Fourier coefficients of 

j)c(/) , representing the envelope functions. According to the discrete time-variable version of 
the Fourier analysis. These harmonic coefficients are given by 27 

■ ' C *k-Tc 

tjc 

Xmc(t)~ \x{t)cQ$a> 0 tdt 
C t k -T c 

2 % 

*™(')— }x(t)smm 0 tdt 

(36) 

where T c is the carrier period and t lz is the time instant. Thus in order to detect the 
amplitudes of the envelope parts of the input signals, it is necessary at first to remove, by 
means of analogue and/or digital initial filtering, the higher harmonic terms and high frequency 
noise from Eq. (36). Then the DFA analysis according to Eq. (37) can be performed. It is easy 
to notice that this algorithm is similar to the synchronous detection method described in 
section 3.1, instead of low-pass filtering after frequency conversion, the DFA algorithm has 
been used. The implementation of this method 27 has shown its very good accuracy and timing 
characteristics. However, a large number of samples in one carrier period is also required 
because of a numerical integration. This method also needs the carrier frequency 
determination. Thus some degree of the carrier frequency shift must be expected and the 
analysis similar to that one presented in section 3.1 ran be performed. 

4. FREQUENCY DEMODULATION ALGORITHMS 
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The frequency demodulation procedure is connected with the computation of an 
instantaneous frequency f(t) of the signal x(t). According to Eq. (1) this operation can be 
written as, 

(37) 

where &(t) is the total instantaneous phase of the signal and /o is the carrier frequency. 
Since this operation requires signal differentiation, in many cases frequency demodulation is 
replaced by phase the demodulation procedure which gives 

J3<f>(fyd{t)^t 
(38) 

where & 0 ^2xrf 0 , Eq. (38) shows that phase demodulation requires removing of the 
constant rotational term & 0 t related to the carrier frequency. Thus the required carrier 
frequency can be determined. Both approaches are similar from the diagnostic point of view 
since they give similar information about the fault. In contrast to amplitude demodulation, 
frequency (phase) demodulation is a nonlinear operation because of the nonlinear 
relationship between the signal x(t) and its instantaneous frequency (phase). Thus the 
procedure is much more difficult. 

There exist many different algorithms of frequency and phase demodulation. Almost 
all of them are inciuded in one of the following categories 26 : (a) FM-to-AM conversion, (b) 
phase-shift discrimination, (c) zero-crossing detection and (d) frequency feedbacks signal 
processing implementation of these algorithms are presented below. 

4.1 . FM-to-AlVI Conversion 



i 
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Any device whose output equals the time derivative of the input produces frequency 
modulation-to-amplitude modulation conversion 26 . Consider the frequency modulated signal 

(39) 



Then 



^ 2.7C at J 



(40) 

Thus taking into account Eq. (38) one can obtain 



2tzA 

(41) 



where Env\x(t)\ is the envelope function of the differentiated signal x(t). Input signal is 

assumed to have some constant amplitude A . In practice an amplitude limiter is necessary at 
the input to remove any variations. Such a limiter is not easy to implement The envelope 
function can be calculated by means of any algorithm presented in section 3. Additionally 
signal differentiation is required. This can be avoided* by analogy to section 3.2.2, by 
measuring acceleration and integration to obtain velocity. This method has been applied in 
machinery diagnostics 10,32 . 



4.2 Synchronous Detection 



Synchronous detection used in section 3.1 to compute amplitude modulation signal 
can be 1 by analogy, used also in the case of phase modulation. Considering 
*(r)=^cos[<3>o**^( f )] and multiplying both sides by 2cos<yo* one obtains 
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2x(/)cos «a?0? = ^ cos f$ (f )+ -4 cos[ 2^ 0 ^ + ^(*) ] 

(42) 

The second term can be, by analogy to section 3.1. removed by low-pass filtering. i 
Thus the instantaneous phase process is easy to extract. The method used simultaneously 
for amplitude and phase demodulation purpose, is very often referred to as complex 

demodulation. Thus^the modulated signal is multiplied by 2e~ J<0 ° t , giving 

(43) 

I 

u Then the procedure to obtain demodulated signals fall onto the following steps 39 : 

Ft (a) determine carrier frequency from the spectrum and bandwidth Ad> occupied by side- 

Si bands, (b) multiply time series by 2e~ JeDot , (c) apply a digital low-pass filter and (d) compute j 

SJ the amplitude and phase of the complex signal obtained from 3. 

=y 

[U More details about this method can be found in 40 . Either in communication systems or 

O In signal processing complex demodulation has received relatively little attention. The main 

lI difficulty is related to the carrier frequency determination. The carrier frequency asynchronous 

destroys harmonic relations in demodulated spectra, 

4.3 Hilbert Transform Approach j 

The analytic signal used in section 3.2.4 to define the envelope function is the basis 
of this approach. Going back to Eq, (22) one can write 

XA (th X (t) + J XH ( t )=a( t )ert<) 

(44) 

i 
i 

i 

i 
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where a(t) is the envelope described in section 2.2.4 and &(t) is the total instantaneous 
phase of the signal being, 

e(t)= arctan + 0 (t) 

(45) 

Taking into account the method of the digital amplitude demodulation described in 
section 3.2.4 k the phase demodulation can be realized. Both operations can be realized 
simultaneously (Fig,4). This algorithm is well established and often used in practice ' ' • . 
The method requires removing the carrier frequency part from the total instantaneous phase. 
This can be done by translating the spectrum components to the negative frequencies, so that 

the component originally at — will be moved to zero. Since in practice the carrier frequency 

is determined with some error instead of #(t) one receives <f>fy+^^-. Using the 

lit ■ 2rt 

Fourier transform and multiplication theorem of the Fourier transform 28 , one will receive in the 

frequency domain the spectrum of the instantaneous phase 



(46) 



which gives an unwanted component in the low frequencies. There exist three methods to 
improve the results. First, the differentiation operation of the total Instantaneous phase d(t) 
can be performed which gives us the instantaneous frequency fit). This operation is not 
effective at all. The other two ways of removing low-frequency trends are 33 least-squares 
polynomial trend removal and high-pass filtering, Since the phase modulation processes are 
low-frequency type, using high-pass filtering, one should be very careful in order not to 
remove the required data. In practice the polynomial trend removal is preferable 33 . 



4.4. Zero-crossing Detection 
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The zero-crossing frequency can be described on the basis of Rice's frequency 
defined as 19 



J/ 2 ^{/)# 



(47) 



From the physical point of view this value indicates the central frequency on which 
the whole power is concentrated. For a narrow-banded Gaussian process Rice's frequency is 
equal to the expected rate of zero-crossings with positive slope, which can be written 

E[nJq)U lim — 

T-±coAT 

(48) 

where N is a number of "positive" zero-crossing values and At is the interval of time. In 
practice Eq. (48) can be replaced by 

*(Ar)=M 
7 At 

(49) 

Thus practically, it resolves into calculations of the interval of time for which the 
change of the phase has a 2k value, which gives us an estimated instantaneous frequency. 
This method used in VA diagnostics 21 requires a great number of samples in one carrier 
period and does not give sufficient accuracy, However, it does not require the carrier 
frequency determination. 



4.5. Wigner-Vilte Distribution 
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The Wigner distribution (WD) can be derived by generalising the relationship between 
the power spectrum and the autocorrelation function for non-stationary, time-variant 
processes. The physical interpretation of the generalised power spectrum Ffaf) is that it 
represents the instantaneous power density spectrum. This leads to the WD defined as*\ 



2 



(50) 

It can be shown that the first derivative of an arbitrary signal's phase reflects the 
mean instantaneous frequency. The WVD can also be used to extract the instantaneous 
frequency: 

Let x(t)-A(t)e 1 ^ where A(t) and #(t), the magnitude and phase respectively are real 
valued functions, the first derivative of the phase is given by 41 , 



LTV J 

(51) 



This simply means that at time t the mean instantaneous frequency of the signal is 
equal to the- mean instantaneous frequency of the WVD. The disadvantage is that at any time 
instant t there is more than one frequency component present, i.e. the instantaneous 
frequency is not a single value, signal energy spreads with respect to the mean instantaneous 
frequencies. 



4.6. Wavelet Transform: Zero-crossings Detection 
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The wavelet transform is used to decompose a signal x(t) into wavelet coefficients 
using the basis of wavelet functions y a ^ (t). This can be expressed as, 

-a: 

(52) 

where is the complex conjugate of ^ (.). There are many functions that can be used 



as the wavelet basis functions, an example is the Morlet wavelet' 



.42 



The square of the modulus of the wavelet transform can be interpreted as an energy 
density over the (a,b) time-scale plane. The energy of a signal is mainly concentrated on the 
time-scale plane around the so called ridge of the wavelet transform. 

It can be shown that 1s,15 r 

\ a J 

(53) 

where & x and ^ denote the instantaneous phases of the signal and the wavelet transform 

respectively. The ridge of the wavelet transform is directly related to the instantaneous 
frequency of the signal The ridge is defined as, 

&(o) 

(54) 

This is used to obtain the instantaneous frequency directly from the ridge of the 
wavelet transform. There exist two algorithms based on the amplitude and phase of the 

IT 
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transform 43 . Application examples include damping estimation procedures" 14-45 and 
identification of nonlinear systems 44 . 

5. EXAMPLES AND APPLICATIONS 

Two Simple examples are shown here to illustrate the application of the instantaneous 
characteristics. The first example shows the power of the Hilbert transform for damage 
detection in gearboxes. The second example from the area of system identification shows 
that often wavelet analysis is better when the instantaneous frequency rather than phase is 

* required in calculations. 

j 5,1 Fault Detection in Gearboxes 

4 

a Fault characteristics in gearboxes appear in the modulated nan-stationary form of 

impacts, which are exhibited in the envelope and instantaneous phase of the vibration data. 
!: The first example involves the analysis of vibration data from a spur gear. The data came 

* from a simple test rig comprising an input gear with 24 teeth driven by an electric motor and 
meshing with 16 teeth of a pinion 13 . The rotational frequencies of the wheel and the pinion are 
25 and 37.5 respectively. This results in the meshing frequency of 600Hz. The 
analysed fault is the loss of part of the tooth due to breakage at a point of the working tip. 
Simply, 1 mm of the facewidth was completely removed. Fig. 2 shows an example of power 
spectra representing normal and damage conditions. The damage condition spectrum 
displays a clear pattern of sidebands around the meshing harmonics. Examples of signal 
demodulation results are given in Fig, 3 and Fig. 4. Time domain averaged signal for damage 
and normal meshing vibration are given in Fig. 3a and Fig. 4a, respectfvely. Local tooth fault 
exhibits an impulse in the envelope function, as shown in Fig. 3b. Also, the change of phase 
can be observed in the instantaneous phase function in Fig. 3c, In contrast, the instantaneous 
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characteristics for the normal meshing vibration remain smooth and do not display any 
disturbances, as shown in Fig. 4b and Fig. 4c. 



5.2. Identification of Nonlinear Systems 



Fig. 5 shows the impulse response function and its spectrum for the vibration 
response from a single degree of freedom rig nonlinear test rig with a cubic stiffness 
characteristic 46 . 

The amplitude of the wavelet transform is given in Fig. 6. The ridge of the wavelet 
transform was computed to obtain the envelope and instantaneous frequency functions for 
the vibration response. Finally, the backbone curve of the analysed system was constructed 
to give the result shown by the solid line in Fig. 7. This characteristic clearly displays cubic 
stiffness nonlinearity. The dashed line in Fig. 7 gives the backbone curve calculated using the 
Hilbert transform approach. Although the oscillations, which are due to differentiation 
procedure, can be removed, the results clearly show the advantage of the wavelet analysis 
over the classical Hilbert transform approach. 



S. CONCLUSIONS 



The envelope and instantaneous frequency/phase functions, are often used for 
damage detection in rotating machinery, for identification of nonlinear systems and in 
deterministic approach to stationarity. 



A number of algorithms for the amplitude and frequency/phase demodulation have 
been briefly described. Finally we are in position to make a comparative analysis. The points 
to be compared are: 

• initial filtering - affects the narrow-band characteristics of the signal; 
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• determination of the carrier frequency - important since the carrier frequency is not 
always known or determined with sufficient accuracy; 

• integration or differentiation - both operations require a sufficiently high sampling 
frequency in order to obtain the required accuracy; 

• initial amplitude limitation - difficult in numerical implementation. 



The comparative resufe are presented in Table 1„ The Hilbert transform and energy 
distribution algorithms seem to be very attractive. The Hilbert transform method is well 
established and easy to implement using the classical FFT algorithm. However, it is only valid 

for single component signals. Also, it requires numerical differentiation for the instantaneous 
frequency, which is not an easy task. New developments in the area of wavelet analysis can 
overtake these drawbacks but often lead to expensive computations. Altogether, it appears 
that the application and implementation very much depend on the problem under 
consideration. 

The envelope and instantaneous frequency/phase functions, are often used for 
damage detection in rotating machinery, for identification of nonlinear systems and in 
deterministic approach to stationary The paper also shows that the envelope and 
instantaneous frequency functions are an important link between classical Fourier approach 
based methods and time-variant analysis. 
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Figure 1„ 

Amplitude and phase modulation signal and its spectrum characteristics. 
Figure 2. 

Power spectra for spur gear vibration data: (a) normal condition (b) damaged tooth. 
Figure 3. 

Signal demodulation characteristics for spur gear vibration data representing normal 
condition: (a) time domain averaged signal (b) envelope (c) instantaneous phase. 

Figure 4. 

Signal demodulation characteristics for spur gear vibration data representing damaged tooth: 
(a) time domain averaged signal (b) envelope (c) instantaneous phase. 

Figure 5. 

Impulse response function (a) and its spectrum far the analysed nonlinear system. 
Figure 6. 

Wavelet transform of the impulse response function. 
Figure 7. 

Backbone curve of the nonlinear system estimated using wavelet analysis («- — ) and the 
Hilbert transform ( — ). 
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Table 1. 



A summary of signal demodulation algorithms. 
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Amplitude [units] 
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Carrier 
frequency 

required ? 


integration 
required? 


Differentiation 
required ? 


Filtering 
required ? 


Comments 


Complex 
Demodulation 
Algorithm 


Yes 


No 


No 


Yes 


« 


Local maxima 
distribution 


No 


Yes 


No 


No 


Problem with time 
of averaging 


Energy 

Distribution 
envelope 


No 


No 


Yes 


No 


Differentiation can 
be replaced by 
integration 


Hilbert 
transform 
algorithm 


No 


No 


No 


No 


Well established 

uses the FFT 
aldorithm but valid 

only for single 
component analysis 


DFA 


Yes 


Yes 


No 


No 




WignerVIHe 
Distribution 


No' 


No 


No 


Yes 


very sensitive to 
noise in the data 


Zero-Crossings 


No 


No 


No 


Yes 


poor resolution and 
accuracy 


Wavelet 
Transform 

j 


No 


No 


No 


No 


suitable for 
multicomponent 

signals but 
computationally 

expensive 
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